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Abstract. Success of precision studies in W and Z boson production at the Tevatron and LHC 
depends on the progress in reducing theoretical uncertainties to the level of 1 — 2%. I review recent 
■ developments in theoretical understanding of W and Z boson production and remaining issues to be 

| solved to meet demands of the near-future experiments. 

o 

Production of W and Z bosons at pp and pp colliders is well-suited for precision 
tests of hadronic matter and electroweak interactions, due to its relatively simple mech- 
^ ; anism and clean experimental signatures. W ± and Z° bosons will be mass-produced in 

the Tevatron Run-2 and at the Large Hadron Collider (LHC) with the goal to measure 
electroweak parameters, probe the internal structure of nucleons, and monitor the col- 
lider luminosity. Deviations from the predictions of the Standard Model may possibly 
occur in the W and Z boson data as a result of new gauge interactions, supersymmetry, 
particle compositeness, or extra spatial dimensions. Searches for new physics are not 
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feasible without adequate evaluation of the large irreducible background of W and Z 



bosons from the conventional processes. 



The measurement of the W boson mass Mw and width Tw is an essential com- 
ponent of the Tevatron physics program. In the Standard Model, radiative corrections 
relate Mw to the square of the top quark mass, mj, and the logarithm of the Higgs boson 
mass, log(M#) [l]. 2 A high-precision measurement of Mw and m t constrains the range 
of Mh allowed within the Standard Model. For example, if Mw and m t are known within 
30 MeV and 2 GeV, respectively, M# is predicted within 35%. Conversely, comparison 
of values of Mw, m t , and Mh measured in different experiments and channels may re- 
veal signs of new physics. The Run-2 plans to measure Mw with accuracy ps 30 MeV 
per channel per experiment, and Tw with accuracy 25-30 MeV [2] . The goal of LHC is 
to reduce 8Mw below 10-15 MeV in order to compete with the precision measurement 
of Mw at the future Linear Collider [3]. 

At collider luminosities above 1 fb _1 , theoretical uncertainties (currently several per- 
cent of the cross section) become the dominant source of error in W and Z observables. 
These uncertainties must be reduced in order to measure the total cross sections with ac- 
curacy < 2% and Mw with accuracy < 0.04%, as envisioned by the Tevatron and LHC 
programs. Several dynamical factors contribute at the percent level, including QCD ra- 
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diative corrections of order 0(af), electroweak corrections of order 0(a), uncertainties 
in parton distributions, and power corrections to resummed differential cross sections. 
Tangible progress has been made in exploring these complex factors, mostly by consid- 
ering each component separately from the rest in order to simplify the problem. Corre- 
lations between effects of different dynamical origins have been often neglected, despite 
their probable significance in the full result. The percent-level QCD and electroweak 
contributions will have to be integrated into a future framework that would account for 
their correlations and be practical enough to serve multifaceted purposes of W and Z 
physics. In this talk, I review the recent advancements and future steps towards develop- 
ment of such a comprehensive model. 

Overview of the process. At hadron colliders, the massive electroweak bosons are 
produced predominantly via qq annihilation and detected by the decay into a pair of 
leptons. A typical W ox Z observable is affected by the QCD and electroweak radiation 
from the hadronic initial state; production, propagation, and decay of massive bosons; 
and electroweak radiation from the leptonic final state. A hadronic cross section is evalu- 
ated in perturbative QCD as a product of hard scattering cross sections and nonperturba- 
tive parton distribution functions (PDF's). The hard scattering cross sections are known 
at &(a%) for inclusive observables (depending on one energy scale), and at <ff{a s ) for 
observables depending on several energy scales. Phenomenological parametrizations of 
the PDF's can be taken from a global fit to hadronic data. Uncertainties in the PDF's 
arising from diverse experimental and theoretical sources can be propagated from the 
global analysis into the predictions for the W and Z cross sections [4]. When the soft 
QCD radiation is enhanced in the scattering (e.g., when the transverse momentum q T 
of the electroweak boson is much smaller than its invariant mass Q), the finite-order 
perturbation theory has to be replaced by resummation of large logarithmic corrections 
through all orders of a s . The resummation is achieved by modifying the form of QCD 
factorization in the affected regions in order to absorb the problematic soft logarithms 
into exponential Sudakov form factors. 

The full 0(a) electroweak corrections are known both for W boson production [5] 
and Z boson production [6]. They depend on the flavor of produced leptons. Not all 
G(a) corrections are of equal importance throughout all phase space. The resonant 
production of W bosons at Q xs Mw is described well within "the pole approximation" 
[7], which neglects the subleading WZ box diagrams. For its part, the pole approximation 
is dominated by virtual corrections to boson propagators and vertices, as well as by QED 
radiation off the final-state charged leptons [8]. The QED radiation from the initial-state 
quarks and interference between the initial and final states are of a smaller magnitude 
numerically, but may have some impact on precision measurements. Away from the W 
boson resonance (Q ^> My/), the pole approximation becomes unreliable and has to be 
replaced by the full 6(a) result [9]. Yet higher-order electroweak corrections may be 
non-negligible as well, especially if enhanced by large logarithms near the boundaries 
of available phase space or at large boson virtualities. W and Z boson production with 
radiation of two additional photons has been studied in Ref. [10]. Leading contributions 
from radiation of multiple soft photons can be resummed through all orders of a, as it 
was done recently for W decays in Refs. [11, 12]. 
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Figure 1. Effect of the &(a%) (NNLO) corrections on the total cross sections a tot (pp — > (W — > lv)X) 

and Ototipp — ► (Z — ► at the Tevatron and LHC [13]. The narrow-width EW approximation was 

used. 



Total cross sections. The interplay of various factors can be most easily demonstrated 
on the example of the total cross sections, given in perturbative QCD by 

<hot{pp ^(V^£i£ 2 )X) = £ / dxidx 2 f a /p(xi,fi)f b j(-){x2 ) fi) 

a,b=q,q,g 

x d tot (ab^{V^hh)X). (1) 

Here V = W or Z; £1^2 are the charged lepton and neutrino (£v) in W boson production, 
or the charged lepton and antilepton (££) in Z boson production; o tot are the hard scatter- 
ing cross sections, known to &(af) [14, 15]; and f a / p (x,n) are the parton distributions. 
The integration is over the partonic momentum fractions x\ and X2, and the summation 
is over the relevant parton flavors. The W and Z total cross sections have been evalu- 
ated recently by the MRST group by using the parton distributions from a recent global 
analysis, realized with inclusion of all known next-to-next-to-leading-order, or NNLO, 
corrections [13]. The &(aj) correction was found to increase a tot by about 4% at the 
Tevatron and reduce o tot by about 2% at the LHC [cf. Fig. 1]. The dependence of o tot 
on the arbitrary factorization scale jU is reduced at 0(oc}) to about 1%, suggesting high 
stability with respect to QCD corrections of order and beyond. 

The magnitude of 0(a) electro weak corrections is comparable to that of the QCD 
corrections. In the limit of the vanishing boson's width Ty, the production of the vector 
bosons can be separated from their decay, by introducing on-shell production cross 

sections o to t(pp — > V) and branching ratios Br(V — > l\£2): 

Ototipp ^(V^ £i£ 2 )X) = o tot (pp ] - V) ■ Br(V - l x £i). (2) 

At this level of accuracy, the branching ratio defined via Eq. (2) coincides with the partial 
width, Br(V — > = r(V — * £\£2)/Ty. The partial width can be defined beyond the 
Born level as a universal pseudo-observable, included in concrete cross sections via 
process-specific relations. For instance, the partial widths for Z bosons are incorporated 
in the fit to the LEP Z-pole observables with the help of a complicated procedure, 
which accounts for the Z line shape, initial- and final- state electroweak radiation, and yZ 




Figure 2. Correlations between the total W and Z cross sections at the Tevatron and LHC found in the 
new CTEQ analysis [17]. The cross sections are computed at ff(a s ) and rescaled by the NNLO /^-factor, 
extracted from Ref. [13]. The dots correspond to 41 PDF sets from the CTEQ6.1 analysis [18]. 



interference (see "The Z boson" review in Ref. [16]). An analogous (but not identical) 
procedure relates the partial widths to the Tevatron cross sections. A tree-level estimate 
in the Standard Model yields T(W -> tv) /T w = 11% and T(Z -> II) /T z = 3.3%, while 
the Review of Particle Physics quotes (10.68 ±0.12)% and (3.3658 ±0.0023)% based 
on the fit to the world (predominantly LEP) data [16]. The shown RPP values are 
averaged over three lepton generations. Consequently the narrow- width approximation 
may deviate from the true Tevatron cross sections by up to ~ 3%, if the Tevatron- specific 
&(a) corrections are not taken into account. 

In addition, the narrow-width approximation neglects the Breit-Wigner line shape 
of vector bosons and is not suitable for computation of differential distributions. The 
next level of refinement is realized in the effective Born approximation (EBA), which 
includes the final width and radiative corrections to the electroweak couplings and 
boson's propagator (but no effects of real particle radiation). With (without) lepton 
identification requirements, the EBA W cross section in Run-2 exceeds the full 0(a) 
cross section by 2% (4%) in the ev e decay channel, and 6% (2%) in the /ivyt decay 
channel [9]. As can be seen from this discussion, the electroweak corrections may 
induce sizable variations in o tot . A percent-level comparison of the total cross sections 
is only meaningful if the employed electroweak corrections are well-documented and 
consistent. 

Uncertainty in o to t due to imprecise knowledge of parton distributions can be eval- 
uated using one of the available methods [4]. The exact values of the PDF errors de- 
pend on the conventions adapted by each group, but all lie in a few-percent range. The 
CTEQ estimate (based on the most conservative criterion among all groups) quotes the 
PDF error of 3% (4%) at the Tevatron (LHC). The PDF errors for W and Z total cross 
sections are highly correlated [cf. Fig. 1], reflecting the striking similarity between the 
quark-dominated initial states of the two processes. The ratio of the two cross sections, 
Rw/z — Otot (W) / Otot (Z), is essentially invariant with respect to the changes in the par- 




ton distributions. A measurement of Z boson cross section can tightly constrain the PDF 
input in W boson production and other kindred processes. Furthermore, fluctuations in 
the Z boson rate in the course of the collider run would reflect changes in the collider 
luminosity, which can be controlled through the measurement of a tot (Z) to the level of 
1 — 2%. These nice features of the W and Z boson total cross sections advance them as 
likely "standard candle" monitors of the partonic and collider luminosities at the LHC 
[19]. Of course, all percent-level effects must be put together in one calculation to realize 
this goal. To reconstruct the total cross section from the visible events, an accurate model 
of the detector acceptance is needed. The W boson acceptances can be controlled at the 
level of 2% by incorporating the parton-level NLO cross sections into parton showering 
programs [20]. 

Rapidity distributions of massive electroweak bosons have been recently computed 
at 0(<xj) with the help of a novel technique based on the optical theorem and recursive 
reduction of arbitrary loop integrals to a small set of known "master integrals" [21]. In a 
large range of rapidities (\y\ < 2), the 0(af) correction leads to essentially uniform 
enhancement of the cross section by 3 — 5% in Z production and 2.5 — 4% in W 
production [cf. Fig. 3]. The corrections are augmented and less uniform in the 

forward rapidity regions. Similarly to the total cross sections, the scale dependence is 
reduced in 0(ocf) rapidity distributions below 1%. 

The shape of the rapidity distributions is determined by the x dependence of parton 
densities. W boson production is sensitive to the flavor composition of the PDF's because 
of the mixing of quark flavors in Wqq coupling. At the pp collider Tevatron, W bosons 
predominantly probe the distributions of u and d quarks (recall that the distribution of 
d antiquaries in an antiproton is equal to the distribution of d quarks in a proton). The 
CTEQ and MRST analyses include the Tevatron W cross sections in the form of the 
charge asymmetry 

, ( do w+ /dy ( -do w ~/dy e 
Ach{yt)- da w +/dye + da w- /dyi ' 

defined in terms of the rapidity distributions do w± /dye of high-p^ charged leptons from 
W ± boson decays [22]. At forward yg, the Tevatron charge asymmetry probes the ratio 




Figure 4. The CDF Run-2 charge asymmetry vs. the resummed result from the qj resummation program 
ResBos [23, 24]. The bands reflect the PDF uncertainty, evaluated using the CTEQ6.1 PDF set [18]. 



d(x, Mw) /u(x,Mw) at* up to about 0.3. In the experiment, selection cuts are imposed on 
the transverse energies Ej of the charged leptons and neutrinos as a means to suppress 
the backgrounds and discriminate between different ranges of x. The preliminary CDF 
Run-2 asymmetry is in a good agreement with the resummed theory cross section at 
25 < E e T < 35 GeV [Fig. 4(a)], but some points fall out of the CTEQ PDF uncertainty 
band at 35 < < 45 GeV [Fig. 4(b)], indicating the need in further improvements of 
the PDF model. 

If the transverse momentum qj of W boson is neglected, a simple kinematic rela- 
tionship exists between the rapidity yi and transverse momentum pn of the charged 
lepton. This relationship induces a strong correlation between the shapes of A c ./,(^) and 
do I dp j so that a measurement of A^iyi) reduces the PDF uncertainty in do /dpji 
and the measured W boson mass [25]. While being a helpful approximation, this corre- 
lation does not hold exactly, due to non-zero qj of W bosons and acceptance constraints 
on the leptonic momenta. If no acceptance cuts are imposed, the charge asymmetry is 
essentially invariant with respect to the QCD corrections [Fig. 5(a)]. The acceptance cuts 
re-introduce the sensitivity of A c ] x (y^) to QCD corrections, due to the differences in phase 
space available at different orders of perturbative calculation. In Fig. 5(b), differences 
can be noticed between the leading-order, next-to-leading-order, and resummed cross 
sections, caused by different shapes of qj distributions in the three calculations. The 
differences are augmented at forward rapidities, where the convergence of perturbation 
series is reduced by enhancement of near-threshold radiation. These features suggest 
that the &(cc^) corrections cannot be replaced by a constant ^-factor at forward rapidi- 
ties. An additional study may be needed to evaluate the impact of QCD radiation and 
phase space restrictions on the observed A^iyi). 

Pt distributions and extraction of A%. The line shape of W bosons is not observed 
in the leptonic channels because of the missing information about the longitudinal mo- 
mentum of the neutrino. For this reason, the W boson mass is commonly deduced from 
the distribution of the transverse momenta of the decay leptons. The distribution in the 

leptonic transverse mass, Mj V = J {\prt\ + |prv|) 2 — (Pre + Ptv) 2 , was the preferred 
observable to extract M w in Run-1 because of its reduced sensitivity to the mechanism 
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Figure 5. (a) Charge asymmetry A c f, (y) of W boson rapidity distributions, da / dy, at various orders of 
a s (no acceptance cuts are imposed) [21]; charge asymmetry A c h(yt) in the Tevatron Run-1 for qj < 30 

GeV, p e T ' v > 25 GeV in the leading-order, next-to-leading order, and resummed calculations [23]. 



of W boson production. The kinematical (Jacobian) peak in the Mj distribution is lo- 
cated exactly at Mj V = M\y at the Born level and smeared by the W boson's width and 
radiative corrections. The shape of the Jacobian peak is not affected, to the first order, 
by the transverse motion of W bosons, caused predominantly by their recoil against 
the soft QCD radiation. For this reason, extraction of Mw from a fit to the M^ dis- 
tribution is less reliant on precise knowledge of the qj spectrum of W bosons and 
PDF's, which is crucial in the other methods [26]. However, the accuracy of the Mj V 
method suffers from the limited precision in the measurement of the neutrino's trans- 
verse momentum (equated in the experiment to the missing transverse energy, $t). An 
alternative technique extracts Mw from the transverse momentum distribution of the 
charged lepton, do /dpxt- This method does not suffer from the complications asso- 
ciated with the reconstruction of , but, as a trade-off, it requires to precisely know 
the qr distribution. The direct measurement of do/dqr involves reconstruction of f.j 
and has low precision. A better prediction for the qj spectrum can be made within the 
theoretical model by applying the resummation methods described below. The W bo- 
son mass can be also extracted from the pr distribution of the neutrinos, or the ratios 
[do/M e T v (W)}/[do/M^ (Z)] and o tot (W)/o tot (Z) of W and Z transverse mass distribu- 
tions and total cross sections [2, 27, 28, 29]. The systematical error tends to be smaller 
in the last two methods because of the cancellation of common uncertainties in the ratio 
of W and Z cross sections. However, the reduced systematical uncertainty is balanced 
in the full result by a larger statistical error, propagated into Mw from a smaller cross 
section for Z boson production. The uncertainties on Mw will be comparable in all meth- 
ods towards the end of Run-2, and several techniques will be probably employed by the 
experimental collaborations as a way to cross check the systematics. 

Factorization at small transverse momenta. The reliability of the finite-order cross 
section is jeopardized at small qj by incomplete cancellation of soft QCD singularities, 
which leaves large logarithms In" (qr / Q) in the hard matrix elements. A re- arrangement 
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Figure 6. The CSS resummed cross sections in Z boson production at the Tevatron. The curves are 
computed in several models for the CSS form factor W(b) at large impact parameters (b > 1 GeV -1 ): 
(a) W(b) at large b is given by extrapolation of its perturbative part from b < 1 GeV -1 (solid); (b) the 
same as (a), multiplied by a Gaussian smearing term e~ 0M (short-dashed); (c) a phenomenological BLN Y 
form, which shows good agreement with the Run- 1 Z data (dot-dashed) [24] ; (d) an updated Ladinsky- 
Yuan form, which shows worse agreement with the Run-1 Z data (long-dashed) [24]. Note that the 
extrapolation model (curves (a) and (b)) must include a Gaussian smearing term e~ s , with g ~ 0.8 Ge V , 
in order to be close to the BLNY form (and, hence, to the data). 

of the perturbation series cures the instability of the theory at q\ <C Q 2 by summing 
the troublesome qj logarithms through all orders of a s into a soft (Sudakov) form 
factor [30]. The validity of such re- arrangement is proved by a factorization theorem 
in the method by Collins, Soper, and Sterman (CSS) [31]. The resummation in vec- 
tor boson production is a special case of a more general problem, and essentially the 
same method applies to hadroproduction in e + e~ scattering [32], and semi-inclusive 
hadroproduction in deep-inelastic scattering [33, 34, 35]. The CSS formalism automat- 
ically preserves the fundamental symmetries (renormalization- and gauge-group invari- 
ance, energy-momentum conservation) and is convenient in practice. The qj resumma- 
tion can be extended to include effects of particle thresholds [36], heavy quark masses 
[37], and hadronic spin [38, 39]. ResBos [23, 24] is a Monte-Carlo integrator program 
that quickly and accurately evaluates the CSS resummed cross sections in Drell-Yan-like 
processes. 

All small-gr logarithms arise in the CSS method from the form factor W(b) in im- 
pact parameter (b) space, composed of the Sudakov exponential and ^-dependent parton 
distribution functions. The resummed qj distribution is obtained by taking the Fourier- 
Bessel transform of W(b) into qj space (realized numerically in ResBos). The alterna- 
tive approaches evaluate the Fourier-Bessel transform of the leading logarithmic towers 
analytically, with the goal to improve transition from the resummed cross section to the 
finite-order cross section at intermediate qj [40, 41]. The integration over all b in the 
Fourier-Bessel transform introduces sensitivity to the nonperturbative QCD dynamics 



at b > 1 GeV . In W and Z boson production, the nonperturbative terms are strongly 
suppressed by the shape of W(b), as well as by the oscillations in the Fourier-Bessel 
transform integrand at large qr- Nonetheless, mild sensitivity to the \axge-b behavior 
remains at qj < 10 GeV, contributing on the top of the overall shape that is tightly 
constrained by the perturbative contributions. The position of the peak in do/dqr can 
move by several hundred MeV depending on the choice of the nonperturbative model 
[cf. Fig. 6]. The range of uncertainty in the W boson's qr distribution is substantially re- 
duced by requiring the nonperturbative model to agree with the measured qj distribution 
in Z boson production. In Run-1, the uncertainty in M\y due to incomplete knowledge of 
the ^rspectrum of W bosons was estimated to amount to 15-20 MeV. 

While significant effort has been put into the study of W(b) at large b [36, 42, 43, 44], 
none of the existing approaches was able so far to adequately describe the observed Z 
boson distribution without introducing free parameters. 3 However, a fit of good quality 
can be made to both Z boson and low-energy Drell-Yan data once a small number of 
free parameters is allowed in the fit to describe the unknown power corrections [24]. 
Variations of these parameters lead to small changes in the qj spectrum of Z bosons and 
large changes in the low-<2 Drell-Yan distributions. The Z boson data taken alone also 
does not constrain the dependence of the nonperturbative terms on Q and parton flavor. 
Consequently a combination of Drell-Yan and Z boson data imposes tighter constraints 
on the nonperturbative parameters than the Z boson data by itself, pretty much as a 
combination of the inclusive data from various experiments places better constraints on 
the PDF's than each individual experiment. 

The feasibility of such a global fit crucially relies on the universality of the nonper- 
turbative contributions. Whether the universality holds is a question under active inves- 
tigation. Universality of the large-& contributions follows from the factorization theorem 
for the resummed cross sections, which was stated in [31] and recently proved in [45]; 
see also the related discussion in Ref. [44]. A high-quality fit to pj data from diverse ex- 
periments would provide, in principle, empirical confirmation of the universality. The fits 
utilizing the b* anzatz [31] - the simplest model of the nonperturbative terms - can in- 
deed achieve a decent % 2 of 170 (130) per 120 data points for the parameter b max equal to 
0.5 (0.8) GeV- 1 . However, the issue needs further investigation because of the substan- 
tial uncertainties in the interpretation of the fits, mostly associated with the low-energy 
Drell-Yan data. The preferred form of the nonperturbative function is strongly correlated 
with the overall normalizations of the Drell-Yan cross sections, which, in their turn, are 
affected by &{a%) corrections (of order 10— 15% [46]), the PDF's, and experimental 
uncertainties. The true uncertainty in the nonperturbative function is linked to the rest of 
the factors, which are not under sufficient control yet in the Drell-Yan process. This is 
particularly true with respect to the PDF parameters, which in principle should be fitted 
together with the nonperturbative function to obviate correlations between the shape of 
do/dqr and parton densities. 



3 It is worth pointing out that the global features of W and Z q T distributions are uniquely determined 
by the perturbative contributions and agree well in the different approaches. The challenging part is to 
describe mild variations at qj below 10 GeV, which cannot be neglected in high-precision measurements. 
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Figure 7. The G[a) electroweak corrections to the leptonic transverse mass distribution, do/dM^, in 
W boson production [7]. No leptonic detection constraints are imposed. 



Electroweak corrections to Mw are dominated by the QED radiation from the final- 
state charged lepton, which results in some loss of the charged lepton's momentum 
to the surrounding cloud of soft and collinear photons. The final-state QED (FQED) 
radiation changes the extracted value of Mw by shifting the Jacobian peak in the My 
distribution in the negative direction [cf. Fig. 7]. In contrast, the initial-state radiation 
and interference terms mostly change the overall normalization of the Jacobian peak 
and have a smaller effect on My/. The size of the electroweak corrections depends on the 
mode of lepton identification. The collinear photons are merged with the electrons in the 
calorimeter towers, but isolated from the muons in the muon detectors. Consequently the 
corrections in the electron channel are larger than in the muon channel in the absence 
of lepton identification requirements, but smaller once the lepton identification is taken 
into the account. The CDF Run- lb analysis estimates the shift in the W boson mass 
due to the FQED radiation by -65 ± 20 MeV and -168 ± 10 MeV in the ev and 
/iV decay channels, respectively [47]. The errors of 20 and 10 MeV correspond to the 
radiative corrections that were ignored in the Run-1 analysis, such as the &{d) initial- 
state radiation and interference, or radiation of two and more photons. These errors need 
to be reduced to push the total 8M\y down to 30 MeV in the Run-2 measurement. 

The correlations between q T spectrum, partonic structure, and electroweak ra- 
diation may prove to be consequential for the upcoming analyses, given that all three 
are the major sources of uncertainties in the extracted W boson mass. In Run-1, these 
interconnections were mostly discarded, so that, for instance, simulations of the QCD 
radiation and electroweak radiation were performed by independent computer programs 
(ResBos [23, 24] and W/ZGrad [6, 7]) and combined ad hoc in the final output. 

Starting from the Run-2, the correlations between the dominant factors cannot be dis- 
missed. The FQED radiation contributes the bulk of the full electroweak correction, and 
it cross-talks with the initial-state QCD radiation through conservation of momentum 
and spin. The combined effect of the G{a) FQED correction and the resummed QCD 
correction was estimated for the Run-2 observables by using a new computer program 
ResBos-A (ResBos with FQED effects) [48]. The FQED and resummed QCD cor- 



rections to the Born-level shape of the Jacobian peak in the Mj distribution were found 
to be approximately (but not completely) independent. The reason is that the My V dis- 
tribution is almost invariant with respect to the transverse momentum of W bosons, so 
that the QCD correction reduces, to the first approximation, to rescaling of the Born- 
level A/|- v distribution by a constant factor. The relationship between FQED and QCD 
corrections is more involved in the leptonic pj distributions, which depend linearly on 
qr ofW bosons. In the pre channel, the combined effect does not factorize into separate 
FQED and QCD corrections to the Born-level cross section. Additional modifications 
will be caused by the FQED correction to the finite part of the resummed cross section 
(Y -piece) and finite resolution of the detector, which were not considered in Ref. [48]. 

The PDF uncertainties were estimated in Run-1 by repeating the analysis for select 
sets of parton densities, which did not cover the full span of allowed variations in the 
PDF parameters. A more systematical estimate can be realized by applying the new 
techniques for the PDF error analysis. The choice of the PDF set affects ^distributions 
directly, by changing the PDF's in the factorized cross section, but also indirectly, by 
modifying the nonperturbative Sudakov function S^pib) in the resummed form factor. 
For a chosen form of SNp(b), the PDF errors can be evaluated within the Hessian matrix 
method, by repeating the computation of qj distributions for an ensemble of sample 
PDF sets. The variations in the resummed qj spectrum for 41 CTEQ6 PDF sets [49] and 
BLNY nonperturbative Sudakov function [24] are shown in Fig. 8. Depending on the 
choice of the PDF set, da / dqj changes by up to ±3% from its value for the central PDF 
set (CTEQ6M). The variations in the PDF's modify both the normalization and shape of 
do/dqr- Although the changes in the shape are relatively weak at qj < 10 GeV, they 
cannot be ignored when M w is extracted from the p T £ distribution. These results do not 
reflect possible correlations between the PDF's and S^p(b) in the global fit to pj data, 
which may be introduced by the dependence of S^p(b) on the normalizations of the low- 
Q Drell-Yan cross sections. In contrast to the PDF's, SNp(b) does appreciably modify the 
shape of do/dqr and may cause larger shifts in the extracted Mw- As discussed earlier, 
the correlation between free parameters in the PDF's and S^p(b) could be explored 
by performing a simultaneous global analysis of the inclusive cross sections and pj- 
dependent data. The first steps towards realization of such a combined fit are being taken 
within the CTEQ collaboration [17]. 

The theoretical assumptions need further scrutiny as well, to guarantee their up-to- 
date level and explore their validity in the new phase space regions accessible in the 
Run-2 and at the LHC. For example, substantial deviations from the present model 
may occur in the resummed cross sections at momentum fractions x of order 10 2 or 
less, typical for forward rapidities in the Run-2, and for all rapidities at the LHC [50]. 
Such deviations may be caused by enhancement of ln(l/x) terms in the perturbative or 
nonperturbative parts of the resummed form factor, which could lead to the broadening 
of q T distributions at small x. An estimate based on the analysis of small-* semi-inclusive 
hadroproduction at HERA suggests that the qj broadening may be visible in a sample 
of forward-rapidity Z bosons collected in the Run-2. If found at the Tevatron, the qj 
broadening will have a profound impact on W and Z production at the LHC. 

Distinctions between the quark flavors were also neglected in the previous resumma- 
tion studies, with the exception of the explicit flavor dependence in the parton densities. 
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Figure 8. The PDF uncertainty in the resummed qj distribution for W boson production in the Tevatron 
Run-2, evaluated using CTEQ6 PDF's [49] and BLNY form of S NP (b) [24]. The lines show the fractional 
differences between do/dqr for 40 sample PDF sets CTEQ6[N] (N = 1, 40), and the central PDF set 
(CTEQ6M). 



This assumption is likely violated at some level, especially for the charm and bottom 
quarks, whose heavy masses suppress very soft radiation. While the heavy-quark sub- 
processes are rare at the Tevatron, they become important at the LHC, where about 17% 
(24%) of W + (W~) bosons will be produced in scattering of charm quarks. The heavy- 
quark masses can be incorporated into the resummation in an extension of the CSS 
formalism to the massive variable- flavor number factorization scheme [37]. We estimate 
within this extension that the heavy-flavor mass terms will have essentially no effect on 
the extracted Mw in the Tevatron Run-2, but may shift Mw by ~ 5 MeV (depending on 
the size of the heavy-flavor nonperturbative terms) at the LHC [51]. 

Conclusions. The theoretical model for W and Z boson production at hadron colliders 
undergoes rapid development. The &((Xj) QCD corrections are becoming available for 
many W and Z observables. The @(cif) corrections show tiny scale dependence and 
can be approximated in several important cases by a uniform rescaling of the &{(X S ) 
cross section by 2.5 — 5%. The 0(aj) correction cannot be replaced by a constant K- 
factor when the lowest-order contribution enters at &{ol s ) (e.g., at large qj or in leptonic 
angular distributions), or near the edges of phase space (e.g., at forward rapidities). The 
finite-order cross sections have to be improved by resummation in kinematical limits 
with enhanced higher-order radiation, notably at qj <C Q, but also at qj 3> Q, Q S> 
My, x — > 0, or x — > 1. Both formal and phenomenological aspects of the resummation 
formalisms require further investigation. 

The &{a) electroweak corrections are also available, and they should be applied 
together with the QCD corrections to reach the required accuracy. Correlations between 
various dynamical factors (PDF's, qj power corrections, electroweak radiation, etc.) 
may be important in the full context and must be systematically explored. Finally, the 
newly found analytical results must be optimally implemented in numerical simulations, 



which should be both fast and precise to adequately serve various practical uses. 
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